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Nomenclature 


Throughout this paper, natural units are used, where h = c = 1. 

E\ab Total energy of the projectile in the lab frame, in units of MeV or GeV. 

Ti a b Kinetic energy of the projectile in the lab frame, in units of MeV or GeV. 

Plab Momentum of the projectile in the lab frame, in units of MeV or GeV. 

m Projectile mass, in units of MeV or GeV. 

y/s Center of momentum total energy, in units of MeV or GeV. 

p Symbol for proton, 

n Symbol for neutron. 

d Symbol for deuteron (bound state of neutron and proton). 

X Generic symbol for unspecified particles. 

7r Symbol for pion. 

7r + Symbol for positively charged pion. 

7T“ Symbol for negatively charged pion. 

7r ± Symbol for charged pions. 

7T° Symbol for neutral pion. 

a Pion production total cross section, in units of mb. 

o-pp-^x Total cross section for the reaction p + p — ► 7 r + X. 


vi 



List of Figures 

1 Parameterized inclusive cross sections for pion production in proton - proton 
collisions given by equations (1) - (3). The green dotted line is a pp-^+x- The red 
solid line is o-pp^ox- The blue long dashed line is cr pp ^-x- The black dotted 
curve is calculated from |(cr pp ^ 7r +x + ^pp-^-x)- Cross sections are set to zero 

for energies below threshold 6 

2 Same as Figure 1, except only the low energy region is shown 6 

3 Same as Figure 1, except only the low energy region is shown and the vertical 

scale is logarithmic 7 

4 Same as Figure 1, except a higher energy region is shown and the vertical and 
horizontal scales are logarithmic. Nonphysical behavior sets in at about 500 GeV. 7 

5 Log - log plot of the 7r + total inclusive cross section parameterization (green, solid 

line) of equation (1) versus experimental data (triangle symbols) of references 
[3,4,6,7,9,14] 8 

6 Same as figure 5, except for use of linear axes 8 

7 Log - log plot of total inclusive cross section parameterization (blue, solid line) of 

equation (2) for the reaction p + p — > n~ + X versus experimental data (triangle 
symbols) of references [3, 4, 7, 9, 14, ?]. Also shown is the cross section data 
(solid squares) from Abdivaliev [10], for the reaction n + p — > n~ + X, compared 

to the parameterization (blue, dashed line) of equation (9) 9 

8 Same as Figure 7, except for use of linear axes and n+p reaction is not shown. . 9 

9 Log - log plot of total inclusive cross section parameterization (red, solid line) of 

equation (3) for the reaction p + p — > 7r° + X versus experimental data (triangle 
symbols) of references [3, 5, 8, 11, 14, 15, 17, 18, 19, 20] 10 

10 Same as Figure 9, except for use of linear axes, and the p + n — > 7r° + X data 

point of Azimov [8] (solid square at 300 GeV) is now shown. The data points for 
p + p — » 7T° + X of Whitmore [3] (solid triangle at 303 GeV) and Azimov [8] (open 
diamond at 300 GeV) have been slightly shifted horizontally so that they do not 
overlap 10 

List of Tables 

1 Reaction thresholds. 7i a b is the kinetic energy threshold, pi a b is the momentum 

threshold, and yfs is the center of momentum total energy at threshold. The mass 
of the deuteron is 1876.125 MeV. Particle C is the produced particle of interest 
and is the first particle listed on the right hand side of each reaction. X are the 
remaining particles. Note that cross sections for n + p are the same as for p + n. 11 

2 Cross section data for pion production in proton - neutron reactions 12 

3 Cross section data for ^ production in proton - proton reactions 13 

4 Cross section data for 7r° production in proton - proton reactions 14 

vii 



Abstract 


Total cross section parameterizations for neutral and charged pion production in 
nucleon - nucleon collisions are compared to experimental data over the projectile 
momentum range from threshold to 300 GeV. Both proton - proton and proton - 
neutron reactions are considered. Overall excellent agreement between parameter- 
izations and experiment is found, except for notable disagreements near threshold. 
In addition, the hypothesis that the neutral pion production cross section can be 
obtained from the average charged pion cross section is checked. The theoretical 
formulas presented in the paper obey this hypothesis for projectile momenta below 
500 GeV. The results presented provide a test of engineering tools used to calculate 
the pion component of space radiation. 


1 Introduction 

Hadron production in nuclear and nucleon collisions can be a significant source of space radi- 
ation, especially for planetary atmospheres. The most important hadron is the pion, because 
it has the lowest threshold for production. While total production cross sections are not used 
in transport equations, it is nevertheless worthwhile to understand them thoroughly because 
integrated differential cross sections must always equal total cross sections. The aim of this re- 
search is to obtain an excellent parameterization of pion production cross sections from nucleon 
- nucleon and nucleus - nucleus collisions. A good understanding of the energy dependence and 
nuclear mass dependence will then provide a further test of the accuracy of differential cross 
sections, when they are integrated. In addition, total cross sections may be used to generate 
differential cross sections. The method for generating these differential cross sections is outlined 
in reference [1]. The idea is to multiply the total cross section by functions that give the shape 
of the differential cross section. One then multiplies by a normalizing factor, so that the re- 
sult equals the total cross section when integrated. This method is useful if the shape of the 
differential cross section is accurately known, because the integrated differential cross section 
is guaranteed to give the correct total cross section. The total cross section parameterizations 
developed previously [2], were compared to only nine experimental data points all above 10 
GeV. The cosmic ray spectrum peaks near the pion threshold at 280 MeV, and therefore it is 
worthwhile to compare to more low energy data. The present work compares parameterizations 
to over thirty experimental data points ranging from threshold to about 300 GeV. Finally, an 
approximation that one often sees in the literature is that the neutral pion cross section can be 
obtained from the average of the charged pion cross sections. This hypothesis will tested. 
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2 Proton - proton reactions 


The total inclusive cross section for pion production in proton - proton collisions has been 
parameterized as [2] 


^pp— >7r+ X 


^pp— >7T _ X 


^pp— >7T°X 


0.00717 + 0.0652 lQ i Tlab + °‘ 162 


-l 


Tisb 

r nnnd „ 0.0846 0.577 N 

0.00456 T n r — T 


7 ^ 0.5 
1 lab 


Hab 


log Tj a K 0 . 3 

0.007+0.1 -g- la - + 


rp2 
1 lab 
-1 


-1 


T\ab 


T 2 

1 lab , 


( 1 ) 

( 2 ) 

( 3 ) 


where Ti a b is the kinetic energy of the projectile, which is given in units of GeV to give a in 
units of mb. These cross sections are plotted in Figures 1 - 4 as a function of pi a b, which is the 
momentum of the projectile. These quantities are related via 


7iab = E isih -m= (4) 

where m is the projectile mass. In equations (1) - (3), the cross sections are set equal to zero 
for energies below threshold. The thresholds are shown in Table 1. 

An interesting aspect of pion production is that the following relation holds. 

^pp— >7r°X ^ 2 ( Cr PP _>>7r-x ^pp— »^r + x) ? (®) 

which is shown as the black dotted curve in Figures 1-4, which closely follows the red solid 
curve (cross section for 7r°), and shows that equation (5) works very well for all energies. Figure 

4 shows that nonphysical behavior sets in at about 500 GeV, where the curves intersect and 
cross over each other. This phenomenon is nonphysical because the cross section should 
not be larger than the 7r + cross section. The parameterizations should not be used above this 
energy. 

The parameterizations in equations (1) - (3) are compared to experimental data in Figures 5 
- 10, which are shown as both linear and logarithmic plots. The reason for this is seen clearly, for 
example, in Figures 5 and 6. Figure 6 seems to show overall very good agreement between the 
parameterization and experiment. However, when plotted on a logarithmic scale, as in Figure 

5 one can now see poor agreement for the data point at 1.38 GeV, which was not apparent on 
the linear plot. 

In reference [2], the parameterizations of equations (1) - (3) were compared only to three data 
points from Whitmore [3] at pi a b = 12, 19 and 24 GeV. In the present work, paramterizations 
are also compared to the Whitmore data at 69, 102, 205 and 303 GeV, and data from other 
authors [4, 5]. More importantly, paramterizations are also compared to some other low energy 
data [6, 7, 8, 9, 10]. Comparison is made to the data point at pi a b = 954 MeV, corresponding 
to Ti a b = 400 MeV, of Abd El-Samad et al. [11]. The cross section for the reaction pp — > pp7r° 
is 0.1 mb. Even though this is an exclusive final state, the exclusive cross section is equal to the 
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inclusive cross section because the energy is below the two pion threshold. See Table 1. The 
complete set of experimental data is listed in Tables 2, 3, and 4. Note that Booth et al. [4] only 
report the cross section for i r + plus ir~ production, defined as a w ± at pi a b = 8.8 GeV. Therefore, 
the individual 7r + and ir~ cross sections listed in Table 3 represent the experimental cross section 
of Booth et al. [4], but with the theoretically calculated individual pion cross section subtracted 
off using equations 1 and 2, i.e. 

< 7 (Table 3) = (Booth) — (equation 2), (6) 

a n - (Table 3) = cr n ± (Booth) — (equation 1). (7) 

Figures 5 and 6 show the comparison of equation (1) to the experimental data. The 
p + p — » 7 r + + X threshold occurs at 7i a b = 289 MeV, or pi a b = 791 MeV, because this is the low- 
est energy reaction for which a 7 r + particle is produced. See Table 1. Comparison between the 
parameterization and experiment is excellent for high energies up to 300 GeV and down to 3.67 
GeV. However, at the lower energy of 1.38 GeV, which is closer to threshold, there is significant 
disagreement between the parameterization and experiment, with with the experimentally ob- 
tained cross sections being much larger than those obtained from the parameterization. Norbury 
et al. [12] have recently shown that cross sections near threshold are important because they 
are enhanced due to larger cosmic ray flux. Figures 7 and 8 show the comparison of equation 
(2) to the experimental data. The p + p — > + X threshold occurs at Ti a b = 600 MeV, or 

P\ab = 1219 MeV (see Table 1). These thresholds are much larger than for 7r + production and 
explains why the 7r + cross section is always larger than the n~ cross section, as shown in Figure 
1. The 7T“ cross section should never be larger than the 7r + cross section because more energy is 
required to produce n~ particles. Comparison between the parameterization and experiment is 
excellent for high energies up to 300 GeV and down to 3.67 GeV. However, at the lower energy 
of 1.38 GeV, there is significant disagreement between the parameterization and experiment, 
this time with the experimentally obtained cross sections being much smaller than those ob- 
tained from the parameterization. Note that now the momentum of 1.38 GeV is very close to 
the 7r“ threshold at 1.219 GeV, and is considered a near threshold cross section [12]. Figures 
9 and 10 show the comparison of equation (3) to the experimental data. The p + p — > 7r° + X 
threshold occurs at Ti a b = 280 MeV, or pi a b = 777 MeV. See Table 1. Comparison between the 
parameterization and experiment is excellent, but there is significant disagreement at the high 
momentum of 300 GeV. The near threshold cross section parameterization is significantly larger 
than experiment. 

3 Neutron - neutron reactions 

Examining the reactions shown in Table 1, reveals that the neutron - neutron (nn) cross section 
for pion production should be the same as for the proton - proton (pp) reaction, if one ignores 
isospin. While Coulomb effects may be important at low energies, the reactions should be equal 
at higher energies. Therefore, 

^nn- ^7rX ~ ^pp— ^7rX (8) 

for all pion species. Unfortunately, there is no data to test this. 
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4 Neutron - proton reactions 

Consideration of the reactions shown in Table 1, plus consideration of multiple pion reactions 
[13] suggests that, at low energy , the neutron - proton (np) cross section might be a simple 
multiple of the pp cross section. The reaction, n + p — > 7r“ + X, has been studied by Abdivaliev 
[10] and the data are shown in Figure 7, compared to the following simple parameterization 
(dashed line) of the np cross section, 

^np— V7T - X 2-5 CTpp^Tj—X 5 (9) 

where cr^^-x is given in equation (2). The np experimental cross sections for other pion 
species (7r + and 7r°) are not available, and parameterizations cannot be developed. 

At higher energies, one may expect pion production cross sections for all pp, np and nn 
reactions to become equal given the proliferation of multiple pion production. This expectation 
seems to be verified by the two data points of Azimov [8] at 300 GeV shown in Figure 10. The 
open triangle is for the reaction p + p — > 7r° + X, and the solid square is for p + n — » 7r° + X. The 
experimental cross section in this high energy region, near 300 GeV, seems to obey the relation 

^pn— ^7T°X ^ ^pp— >7T°X ■ (1®) 

Finally, one would expect 

^np— >-7rX = ^pn— »7rX (H) 

to hold at all energies for all pion species. 


5 Conclusions 

Total cross section parameterizations for neutral and charged pion production in nucleon - 
nucleon collisions have been compared to an extensive set of experimental data over the projectile 
momentum range from threshold to 300 GeV. For ^ production, the disagreement between the 
parameterization and experiment near threshold at pi a b = 1.38 GeV needs to be resolved. This 
disagreement is not obvious in the linear plots, but shows up clearly in the logarithmic plots. 
For 7T° production, the disagreement between the parameterization and experiment at pi a b = 
300 GeV needs to be resolved. There is a need for more experimental measurements of pn cross 
sections for all pion species at low and high energy. In spite of these disagreements, there is 
excellent overall agreement between parameterizations and experiment. Above the energy region 
near threshold, the parameterizations provide accurate cross section predictions. The results can 
be used to validate engineering tools used to calculate the pion component of space radiation. 
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o (mb) 


p + p -> 7T + X 



Figure 1: Parameterized inclusive cross sections for pion production in proton - proton 
collisions given by equations (1) - (3). The green dotted line is CTpp-^+x- The red solid line 
is CTpp-^ox- The blue long dashed line is CTpp-^-x- The black dotted curve is calculated 
from | (crpp-^+x + 0 pp-» 7 r-x)- Cross sections are set to zero for energies below threshold. 
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Figure 2: Same as Figure 1, except only the low energy region is shown. 
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Figure 3: Same as Figure 1, except only the low energy region is shown and the vertical 
scale is logarithmic. 
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Figure 4: Same as Figure 1, except a higher energy region is shown and the vertical and 
horizontal scales are logarithmic. Nonphysical behavior sets in at about 500 GeV. 
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Figure 5: Log - log plot of the 7r + total inclusive cross section parameterization (green 
solid line) of equation (1) versus experimental data (triangle symbols) of references [3, 4 
6, 7, 9, 14]. 
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Figure 6: Same as figure 5, except for use of linear axes. 
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(p,n) + p -> 7r + X 
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Figure 7: Log - log plot of total inclusive cross section parameterization (blue, solid line) of 
equation (2) for the reaction p + p —*• tt ~ + X versus experimental data (triangle symbols) 
of references [3, 4, 7, 9, 14, ?]. Also shown is the cross section data (solid squares) from 
Abdivaliev [10], for the reaction n + p — ► 7r - + X, compared to the parameterization (blue, 
dashed line) of equation (9). 



140 
120 
100 

cr(mb) 80 

60 

40 
20 
0 

0 50 100 150 200 250 300 

piab(GeV) 

Figure 8: Same as Figure 7, except for use of linear axes and n+p reaction is not shown. 
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P + P -> 7T° + X 
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Figure 9: Log - log plot of total inclusive cross section parameterization (red, solid line) of 
equation (3) for the reaction p + p — > 7 T° + X versus experimental data (triangle symbols) 
of references [3, 5, 8, 11, 14, 15, 17, 18, 19, 20]. 
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Figure 10: Same as Figure 9, except for use of linear axes, and the p + n — > 7 T° + X 
data point of Azimov [8] (solid square at 300 GeV) is now shown. The data points for 
p+p — > 7 r°+X of Whitmore [3] (solid triangle at 303 GeV) and Azimov [8] (open diamond 
at 300 GeV) have been slightly shifted horizontally so that they do not overlap. 
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Table 1: Reaction thresholds. 7] a b is the kinetic energy threshold, pi a b is the momentum 
threshold, and y/s is the center of momentum total energy at threshold. The mass of 
the deuteron is 1876.125 MeV. Particle C is the produced particle of interest and is the 
first particle listed on the right hand side of each reaction. X are the remaining particles. 
Note that cross sections for n + p are the same as for p + n. 


A + B — > C + X 7] a b Pi a b y/s 

(MeV) (MeV) (MeV) 


P + P— > 7T° + p + p 

7T° T 7T° T P T P 
7T + + p + n 
7T + + d 

7T“ + 7T + + p + p 

n + n— > 7T° + n + n 

7T + + 7T“ + n + n 
7T“ + p + n 
7T~ + d 

p + n— > 7r° + p + n 
7T° + d 

7r + + n + n 
7T — T p T p 


279.7 

776.6 

2011.5 

578.7 

1192.1 

2146.5 

292.3 

796.2 

2017.4 

288.6 

790.5 

2015.7 

599.8 

1218.7 

2155.7 

279.7 

777.0 

2014.1 

599.8 

1219.3 

2158.3 

286.7 

788.0 

2017.4 

283.0 

782.3 

2015.7 

279.5 

776.2 

2012.8 

275.8 

770.5 

2011.1 

292.1 

795.9 

2018.7 

286.5 

787.3 

2016.1 
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Table 2: Cross section data for pion production in proton - neutron reactions. 


Reaction 

T lab 
(GeV) 

P\ab 

(GeV) 

V? 

(GeV) 

a 

(mb) 

Reference 

p + n — > 7T° + X 

299.1 

300 

23.8 

136 ± 10 

[8] 

n + p — > 7T“ + X 

0.62 

1.25 

2.17 

1.45 ± 0.06 

[10] 


1.03 

1.73 

2.34 

4.1 ± 0.1 

[10] 


1.48 

2.23 

2.51 

8.7 ± 0.3 

[10] 


3.00 

3.83 

3.03 

16.44 ± 0.51 

[10] 


4.25 

5.1 

3.39 

21.0 ± 0.53 

[10] 
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Table 3: Cross section data for ^ production in proton - proton reactions. 


Reaction 

Tlab 

(GeV) 

Plab 

(GeV) 

V* 

(GeV) 

o 

(mb) 

Reference 

P + P — > 7T + + X 

0.59 

1.2 

2.2 

9.7 ± 1.2 

[6] 


0.73 

1.38 

2.21 

13.5 ± 0.73 

[9] 


2.85 

3.67 

2.98 

26.1 ± 5.0 

[7] 


7.91 

8.8 

4.3 

39.9 ± 0.6 

[4] 


11.1 

12 

4.9 

43.2 ± 0.9 

[3, 14] 


18.1 

19 

6.1 

47.5 ± 1.0 

[3] 


23.1 

24 

6.8 

56.8 ± 0.9 

[3, 14] 


68.1 

69 

11.5 

81.5 ± 8.2 

[3] 


101.1 

102 

13.9 

91.8 ± 9.2 

[3] 


204.1 

205 

19.7 

108 ± 11 

[3] 


302.1 

303 

23.9 

125 ± 13 

[3] 

P + P — ► 7T“ + X 

0.73 

1.38 

2.2 

0.03 ± 0.01 

[9] 


2.85 

3.67 

2.98 

6.0 ± 1.2 

[7] 


7.91 

8.8 

4.3 

19.1 ± 0.6 

[4] 


11.1 

12 

4.9 

21.3 ± 0.4 

[3, 14] 


18.1 

19 

6.1 

30 ± 0.6 

[3, ?] 


23.1 

24 

6.8 

33.8 ± 0.6 

[3, 14] 


68.1 

69 

11.5 

62.7 ± 1.0 

[3] 


101.1 

102 

13.9 

66.9 ± 1.3 

[3] 


204.1 

205 

19.7 

86 ± 2 

[3] 


302.1 

303 

23.9 

99.5 ± 3.0 

[3] 
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Table 4: Cross section data for 7 r° production in proton - proton reactions. 


Reaction 

Tlab 

(GeV) 

Plab 

(GeV) 

(GeV) 

a 

(mb) 

Reference 

P + p — > 7T° + X 

0.4 

0.954 

2.07 

0.1 ± 0.007 

[11] 


11.1 

12 

4.9 

31.5 ± 2.6 

[3, 14, 15] 


18.1 

19 

6.1 

41.6 ± 3.5 

[3] 


23.1 

24 

6.8 

53.5 ± 3.1 

[3, 14] 


68.1 

69 

11.5 

81.1 ± 4 

[3, 16, 17] 


101.1 

102 

13.9 

85 ± 8 

[3, 18] 


204.1 

205 

19.7 

107.3 ± 6.8 

[3, 19] 


299.1 

300 

23.8 

130 ± 5 

[8] 


302.1 

303 

23.9 

127 ± 12 

[3, 20] 


399.1 

400 

27.4 

142.3 ± 6.2 

[5] 
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